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Numerical  models  were  developed  to  assess  the  thermal  performance  of  a  HEX  reactor  with  solid-to-gas 
reacting  flow.  Based  on  the  experimental  results  obtained  in  part  I  of  this  study,  numerical  models  were 
developed  to  predict  the  thermal  performance  of  a  plate  heat  exchanger-based  HEX  reactor  involving 
multi-phase  flow  with  chemical  reactions.  A  reduced-order  numerical  model  of  a  chevron  plate  heat 
exchanger  was  developed  with  thermal  and  momentum  transfer  analogies.  Empirical  correlations  for 
momentum  transfer  and  void  fraction  (validated  in  part  I  of  this  study)  were  implemented  in  the  numer¬ 
ical  model.  The  numerical  model,  coded  in  Maple  13™,  was  used  to  size  a  compact  reactor  with  a  thermal 
load  rating  of  2  kW  for  the  desired  operating  temperature;  the  modeling  framework  developed  can  also 
be  used  to  study  different  candidate  gas-generating  reacting  species,  working  fluids,  and  PHE 
configurations. 
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1.  Introduction 

Chevron  PHEs  are  known  to  provide  high  heat  transfer  coeffi¬ 
cients  with  a  high  surface  area-to-volume  ratio  [1],  enabling  high 
thermal  power  ratings  without  incurring  significant  footprint  ex¬ 
penses.  The  corrugated  plate  geometry  also  produces  swirling  sec¬ 
ondary  flows  and  early  transition  to  turbulence  at  low  Re,  yielding 
near-ideal  mixing  behavior  in  addition  to  exceptional  thermal  per¬ 
formance:  for  instance,  Santacesaria  et  al.  [2]  reported  higher 
yields  of  biodiesel  at  relatively  low  Re,  which  they  credited  to 
strong  localized  micromixing  and  turbulence.  Edge  et  al.  [3]  used 
a  model  single-phase  exothermic  reaction  to  characterize  both 
heat  transfer  and  reactant  micromixing  performance  for  a  PHE- 
based  HEX  reactor.  Their  results  showed  reduced  byproduct  forma¬ 
tion  indicative  of  high  micromixing. 

The  results  reported  in  the  literature  suggest  tremendous 
potential  for  the  HEX  reactor  scheme  introduced  in  a  previous 
study  [4],  and  discussed  in  part  1  of  this  study.  Nevertheless,  it  is 
clear  that  there  is  a  lack  of  both  experimental  and  theoretical 
understanding  of  solid-to-gas  reacting  flows  in  a  chevron  PHE. 
Consideration  of  such  a  unique  HEX  reactor  will  require  extensive 
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exploration  from  an  analysis  and  design  standpoint.  It  is  therefore 
important  to  develop  reliable  and  accurate  models  to  predict  both 
the  thermal  performance  and  conversion  efficiency  of  a  given  HEX 
reactor  design,  in  effect  providing  an  initial  feasibility  test  for  the 
HEX  reactor’s  role  as  a  thermal  management  scheme.  There  ap¬ 
pears  to  be  only  a  limited  number  of  models  developed  for  chevron 
PHE-based  HEX  reactors  (such  as  that  developed  by  Edge  et  al.  [3]), 
none  of  which  consider  multi-phase  flows.  On  the  other  hand, 
models  developed  for  chevron  PHEs  that  consider  multi-phase 
flows  focus  on  refrigerant  evaporation  and/or  condensation  (e.g., 
[5,6]),  but  without  chemical  reaction.  To  our  knowledge,  there 
are  no  models  proposed  in  the  open  literature  that  are  applicable 
to  solid-to-gas  reacting  flow  in  a  chevron  PHE.  Hence,  the  focus 
of  this  paper  is  to  develop  a  computer  model  for  the  coupled  heat, 
momentum,  mass  transfer,  and  chemical  kinetics  of  an  endother¬ 
mic  solid-gas  reaction  in  a  chevron  PHE.  The  model  predictions 
for  thermal-hydraulic  performance  were  used  to  demonstrate 
the  utility  of  the  solid-gas  endothermic  HEX  reactor  for  thermal 
management.  Calculations  were  performed  using  the  void  fraction 
and  two-phase  multiplier  correlations  that  were  experimentally 
validated  in  part  1.  Thermal  performance  estimates  were  obtained 
using  a  modified  heat/momentum  transfer  analogy  from  the  liter¬ 
ature  that  has  been  developed  for  PHEs. 

The  volumetric  energy  densities  of  thermochemical  reactions 
generally  exceed  that  of  thermophysical  reactions.  It  has  been 
previously  demonstrated  [4]  that  the  endothermic  solid-to-gas 
decomposition  of  ammonium  carbamate  (AC)  holds  significant 
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Nomenclature 

A 

surface  area 

q" 

wall  heat  flux 

C 

mass  concentration 

Opiate 

plate  thickness 

D 

diameter 

U 

velocity 

Da 

Damkohler  number 

W 

power  rating 

G 

constant  in  Eq.  (3) 

X 

mass  quality 

L 

length 

z 

dimensionless  axial  coordinate 

M'" 

dimensionless  reaction  rate 

z* 

axial  coordinate 

Nu 

Nusselt  number 

P 

pressure 

Greeks 

Pp 

partial  pressure 

A 

corrugation  pitch 

Pr 

Prandtl  number 

© 

dimensionless  temperature  in  Figs.  9  and  10 

Qlatent 

heat  rejected  to  ammonium  carbamate 

area  enlargement  factor 

Qload 

total  heat  rejected 

Ojp 

two-phase  multiplier 

R 

heat  capacity  ratio 

a 

void  fraction 

Re 

Reynolds  number 

P 

chevron  angle 

T 

temperature 

9 

dimensionless  temperature 

r 

normalized  load  temperature 

s 

Qlatent/ Qload 

Th 

load  temperature  at  full  conversion 

K 

rate  constant 

U 

superficial  velocity 

n 

dimensionless  concentration 

V 

volume 

p 

density 

w 

plate  width  (between  gasket) 

X 

Lockhart-Martinelli  parameter 

Subscripts 

A 

acceleration 

Lowercases 

F 

friction 

a 

rate  constant  in  Eq.  (2) 

G 

gravity 

d 

1/2  corrugation  depth 

SP 

single-phase 

Cp 

specific  heat  capacity 

TP 

two-phase 

e 

constant  in  Eq.  (3) 

c 

corrugated  section 

f 

Darcy  friction  factor 

e 

equivalent 

g 

gravitational  acceleration 

eq 

at  equilibrium 

h 

convection  coefficient 

g 

gas  phase 

yheff 

overall  thermal  resistance 

h 

hydraulic 

i 

cell  index 

in 

value  at  inlet 

k 

thermal  conductivity 

Id 

liquid  phase,  load  side 

^^plate 

thermal  conductivity  of  plate 

Is 

liquid  phase,  slurry  side 

fn'" 

volumetric  reaction  rate 

P 

AC  particle 

Ktal 

total  mass  flux 

PP 

port-to-port 

n 

channel  index 

potential  as  a  high  energy  density  thermochemical  heat  sink  for 
thermal  management.  AC  is  a  salt  that  is  readily  available  in  pow¬ 
der  form  as  a  cheap  byproduct  of  urea  synthesis  [7].  AC  undergoes 
a  reversible  solid-gas  decomposition  reaction,  shown  below: 

NH2COONH4(s)  ^  C02(g)  +  2NH3(g)  (1) 

Despite  the  low  operating  temperature  (~70  °C)  and  high  en¬ 
ergy  density  (2000  kj/kg),  the  impediment  to  AC  in  a  practical  ther¬ 
mal  management  system  is  its  low  thermal  conductivity  (~0.4  W/ 
m  K)  [7].  This  presents  a  significant  barrier  to  conducting  thermal 
energy  in  bulk  AC  salt.  Schmidt  [7]  proposed  forming  reacting  slur¬ 
ry  by  immersing  small  AC  salt  particles  in  a  non-reactive  heat 
transfer  fluid  (HTF)  to  overcome  this  impediment.  Previous  re¬ 
search  [8,9]  demonstrated  successful  use  of  PHEs  with  slurries  con¬ 
taining  thermophysical  phase-change  materials  (PCM).  However, 
due  to  the  high  enthalpy  of  decomposition,  AC  can  enable  higher 
overall  energy  density  to  be  achieved  than  with  the  thermophysi¬ 
cal  PCM  slurries  considered  in  the  aforementioned  studies. 

2.  Mathematical  analysis 

Gasketed  chevron  PHEs  are  a  stack  of  channels  formed  by  thin 
corrugated  heat  transfer  plates  which  are  separated  by  gaskets. 
The  corrugations  are  aligned  in  a  repeating  series  of  chevrons  to 


form  an  oblique,  wavy  surface.  A  discussion  of  the  background  of 
PHEs  is  provided  by  Johnson  et  al.  [4]  and  in  part  1  of  this  study. 
The  geometric  characteristics  of  the  PHE  used  in  this  study  are  de¬ 
picted  in  Eig.  1. 

A  schematic  of  the  HEX  reactor  architecture  considered  in  this 
study  is  shown  in  Eig.  2.  A  hot  stream  fluid  delivers  heat  to  the  cold 
stream  fluid  carrying  suspended  AC;  the  PHE  has  a  single-pass, 
parallel,  vertical  upflow  arrangement.  The  model  developed  in  this 
study  considers  the  PHE  to  comprise  of  a  stack  of  equivalent  flat 
plates,  forming  a  collection  of  N  channels  as  illustrated  in  Eig.  2. 
Each  channel  is  represented  by  a  ‘submodel’  whose  interaction 
with  other  submodels  in  the  system  is  defined  by  appropriate 
boundary  conditions,  as  illustrated  in  Eig.  3.  In  the  case  considered 
here,  channel  submodels  are  permitted  to  interact  by  exchanging 
thermal  energy.  The  heat  transfer  coefficient,  friction  factor,  two- 
phase  multiplier,  and  void  fraction  are  estimated  by  employing 
empirical  correlations  obtained  from  the  literature. 

This  approach  is  based  on  the  idea  of  hierarchical  modeling, 
which  has  been  previously  employed  for  the  analysis  of  PHE 
performance  in  industrial  applications  [10].  The  treatment  of 
individual  channels  as  separate  submodels  (Eig.  3)  also  allows  for 
non-idealities  seen  in  real  heat  exchangers,  such  as  maldistribution 
of  flow  patterns,  to  be  incorporated  into  the  overall  reactor  model 
(Fig.  2). 
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Fig.  1.  (Left)  Image  of  a  representative  chevron  PHE  heat  transfer  plate  (Alpha  Laval  MF3)  used  in  this  study.  (Right)  Wavy  corrugation  profile  observed  in  side  view. 


2A.  Transport  equations 


Due  to  the  small  particle  size,  the  contribution  of  the  suspended 
solids  on  the  momentum  equation  is  negligible  [12]  (i.e.,  the 
small  concentration  of  the  tiny  suspended  particles  does  not 
cause  appreciable  change  in  the  effective  viscosity  and  rheolog¬ 
ical  properties  of  the  bulk  liquid  phase). 

Reduced  order  models  involving  one-dimensional  equivalent 
flows  are  used;  i.e.,  cross-sectional  variation  in  the  temperature, 
velocity,  and  void  fraction  profiles  are  replaced  by  area  aver¬ 
aged  quantities. 

0  Primary  flow  variables  are  permitted  to  differ  between  the 
liquid  and  gas  phases  (separated  flow  model). 

0  Heat  transfer  to  the  gas  phase  is  negligible  when  compared 
to  heat  transferred  to  the  liquid  phase  and  the  suspended 
AC.  This  is  a  reasonable  engineering  approximation  due  to 
the  short  residence  times  and  the  low  thermal  conductivity 
values  of  the  gas  phase  compared  to  that  of  the  liquid  phase. 
The  reaction  occurs  at  the  local  bulk  temperature  due  to  homo¬ 
geneous  mixing  produced  by  bubble  agitation  and  secondary 
flows  characteristic  of  PHEs. 

Negligible  conduction  heat  transfer  in  the  axial  direction. 
Fluids  are  Newtonian  and  incompressible. 


The  separated  flow  model  [11  ]  is  used  for  developing  the  trans¬ 
port  equations  of  energy,  mass,  and  momentum  for  the  liquid  and 
gas  phases  in  the  multi-phase  flow  configuration,  wherein  the  gas 
phase  and  liquid  phase  (heat  transfer  fluid  and  suspended  solid) 
are  considered  segregated.  The  mutual  interactions  between  gas 
and  liquid  phase  are  estimated  using  appropriate  experimental 
correlations  (validated  from  the  experiments  conducted  in  part  1 
of  this  study). 

Some  of  the  assumptions  used  in  the  HEX  reactor  model  are 
listed  as  follows: 

•  Steady-state  operation  (using  time  averaged  values). 

•  The  gas  mixtures  (reaction  products)  are  insoluble  in  the  liquid 
phase. 

•  No  evaporation  occurs  in  the  liquid  phase  on  the  slurry  side. 

•  The  flow  is  approximated  as  a  pseudo  two-phase  flow  (homoge¬ 
neous  flow  model  for  the  AC  in  the  liquid)  where  the  liquid 
“solution”  of  AC  (rather  than  a  collection  of  discrete  particles) 
is  considered  to  undergo  uniform  volumetric  heat  generation 
(or  heat  sink)  due  to  homogeneously  distributed  volumetric 
chemical  reactions  and  mass  concentration  at  a  given  axial 
location. 


In  accordance  with  the  aforesaid  assumptions,  an  integral  mass, 
species,  energy,  and  momentum  balance  on  a  1 -dimensional  control 
volume  yields  the  following  governing  equations: 


Plate  j-1 

U(n)  - - 

TJn)  _ V 

CJn) 

PJn)  ^ 


Plate  j 


U(n+1) 

TJn+1) 


Plate  j+1 


1/  \7~  \t 


^  Cold  (Slurry)  Side 
)  Channel  w 


//  .  t.  Hot  (Load)  Side 

^ w,n+\  )  Channel  n+1 


Fig.  3.  Interaction  between  channel  submodels. 
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2.1.1.  Continuity 


2.2.  Incorporation  of  kinetic  model 


Uls,n  = 


Uls,n 

1  -a„ 


(2) 


tlcrri  — 


klls,n(^in,n  ~  ^n) 

Pg,n^n 


(3) 


where  u  is  the  velocity  (m/s),  U  is  the  superficial  velocity  (m/s),  p  is 
the  density  (kg/m^),  and  C  is  the  mass  concentration  of  solid  AC  (kg/ 
m^)  per  unit  volume  of  liquid  slurry. 


2.1.2.  Species 


=-<(!-«„)  (4) 

where  fn'"  is  the  volumetric  reaction  rate  (kg/m^  s)  per  unit  volume 
of  liquid  slurry. 

2.1.3.  Energy 


Uls,n 


din 

dz* 


AH 

Pls^p,ls 


<(1 


Otn)  + 


qn+qn-1 

^ePist^p,ls 


(5) 


U 


dTn 


Q'n  +  q'^1 

DePidCpJd 


(6) 


where  T  is  the  temperature  (K),  q"  is  the  wall  heat  flux  (W/m^),  Cp  is 
the  specific  heat  capacity  (J/kg  K),  and  is  the  equivalent  diameter 
(m,  twice  the  corrugation  amplitude).  The  momentum  equation 
splits  the  total  pressure  gradient  into  frictional,  hydrostatic,  and 
acceleration  components  [11]: 


2.1.4.  Momentum 


dPn  ^  dPF,n  dPG,n  dPA,n 

dz*  dz*  dz*  dz* 

where. 


dPA,n 

dz* 


tlls,nPls(^  ^n) 


duisp 

dz* 


+  ^g,nPgpOCn 


dUgn 

dz* 


(7) 

(8) 


The  current  state  of  knowledge  about  the  kinetics  of  AC  decom¬ 
position  is  quite  limited;  the  open  literature  offers  few  mathemat¬ 
ical  models  for  the  reaction  rate.  As  a  starting  point,  we  selected 
Claudel  and  Boulamri’s  [13]  net  rate  model  for  the  decomposition 
or  formation  of  AC,  as  provided  below: 

dP 

-^  =  aA,\P,-Pe/  (12) 

where  Pp  is  the  partial  pressure  (kPa)  of  the  NH3  and  CO2  evolved 
during  the  course  of  decomposition,  a  is  the  rate  coefficient  (kPa/ 
m^  s),  and  Ap  is  the  surface  area  of  solid  AC  (m^).  The  equilibrium 
pressure.  Peg,  follows  an  Arrhenius-type  dependence  on  tempera¬ 
ture  [14]: 

Pe,(T)  =  Gexp(-|)  (13) 

where  G  and  e  are  constants  equal  to  1.9122  x  10^°  kPa  and 
6321.7  K,  respectively,  as  obtained  from  Koutinas  et  al.  [14].  Fig.  4 
shows  a  plot  of  Eq.  (13)  as  a  function  of  temperature,  with  values 
measured  from  various  sources  in  the  literature.  According  to  Eq. 
(13),  the  minimum  operating  temperature  (under  standard  atmo¬ 
spheric  pressure)  for  the  system  is  ~60  °C.  Schmidt  [7]  also  mea¬ 
sured  the  decomposition  pressure  as  a  function  of  temperature  for 
AC  immersed  in  propylene  glycol  or  ethylene  glycol,  and  found  no 
significant  departure  from  the  curve  described  by  Eq.  (13). 

The  red  arrow  in  Eig.  4  represents  the  thermodynamic  “driving 
force”  behind  the  reaction,  or  “equilibrium  drop,”  as  some  authors 
have  used  for  solid-to-gas  reversible  reactions  [18].  The  red  dotted 
line  represents  the  equilibrium  vapor  pressure  corresponding  to 
the  local  temperature  of  the  solid  AC,  whereas  the  black  dotted  line 
represents  the  local  partial  pressure  acting  on  the  reactor  contents. 
Schmidt  confirmed  that  the  decomposition  rate  could  be  acceler¬ 
ated  appreciably  by  maintaining  the  AC  far  from  its  equilibrium 
condition,  maximizing  the  equilibrium  drop.  This  driving  force  is 
accounted  for  in  Claudel  and  Boulamri’s  model.  Since  the  form  of 
Eq.  (12)  is  indicative  of  a  physical  mechanism  [13],  we  postulate 
that  the  same  mechanism  will  exist  even  when  immersed  in  a  heat 
transfer  fluid. 

The  rate  equation  discussed  above  is  incorporated  into  the  spe¬ 
cies  conservation  equation  via  the  volumetric  reaction  term,  m'". 


dPG,n 

dz* 


g[0  -  a.n)Pis  +  anPg,„l 


(9) 


dPf.n  _(fpU^)ls.n^2 
df  ~  2Du 


(10) 


where  g  is  the  gravitational  acceleration  (m/s^),  is  the  hydraulic 
diameter  (m).  The  subscript  n  is  used  to  indicate  that  the  respective 
quantity  may  vary  from  channel  to  channel.  Inlet  temperatures, 
pressures,  concentrations,  and  flow  rates  to  each  channel  are  pre¬ 
sumed  known.  Heat  transfer  between  channel  submodels  n  and 
n  -  1  is  assumed  to  occur  by  1-D  heat  conduction,  as  is  typical  in 
heat  exchanger  analysis  [1]: 


9n-l  — 


^  _l_  Opiate  ^ 
hn  kpiate  ^n-1 


(Tn-t-Tn)=heff,n-t(Tn-t-Tn)  (H) 


where  h  is  the  local  convection  coefficient  (W/m^  1<),  tpiate  is  the  heat 
transfer  plate  thickness  (m),  and  kpiate  is  the  thermal  conductivity  of 
the  plate  (W/m  K).  The  end  plates  are  treated  as  insulated  bound¬ 
aries,  since  they  are  much  thicker  than  the  plates,  and  PHEs  are  de¬ 
signed  to  minimize  heat  losses  to  the  environment  through  the  end 
plates  and  gaskets  [1]. 


Fig.  4.  Dissociation  pressure  Peg  as  a  function  of  temperature  predicted  by  Eq.  (3) 
compared  with  experimental  results  obtained  in  the  literature.  Red  squares: 
Koutinas  et  al.  [14];  green  triangles:  Bennet  et  al.  [15];  purple  X’s:  Egan  and  Potts 
[16]:  blue  asterisks:  Briggs  and  Migrdichian  [17].  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version 
of  this  article.) 
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Assuming  that  the  dependence  upon  surface  area  and  equilibrium 
drop  is  the  same  as  that  which  was  proposed  by  Claudel  and 
Boulamri  [13],  the  decomposition  rate  of  an  individual  particle  is: 

^^(p^V,)  =  -KA,\Pe,-P\^  (14) 

where  Vp  is  the  particle  volume  (m^),  pp  is  the  particle  density  (kg/ 
m^)  (assumed  to  be  uniform  throughout  the  particle),  and  k  is  the 
specific  rate  coefficient  based  on  surface  area,  which  will  assuredly 
differ  from  the  rate  constant  calculated  in  Claudel  and  Boulamri’s 
study.  As  an  initial  approximation,  since  the  variation  of  NH3  and 
CO2  partial  pressures  in  the  HEX  reactor  has  not  yet  been  experi¬ 
mentally  investigated,  the  partial  pressure  Pp  in  Eq.  (14)  is  assumed 
to  equal  the  local  static  pressure  of  the  liquid,  P.  Peq  is  a  function  of 
temperature;  Schmidt’s  [7]  results  for  AC  in  propylene  glycol  and 
ethylene  glycol  match  well  with  the  Peq(T)  curves  obtained  in  liter¬ 
ature  for  dry  AC  [14],  as  described  by  Eq.  (13). 

It  follows  that  the  mass  concentration  of  solid  AC  per  unit  volume 
of  liquid  heat  transfer  fluid  is  C  =  n"' PpVp.  On  a  volumetric  basis,  the 
rate  of  solid  mass  consumption  by  decomposition  may  then  be  writ¬ 
ten,  assuming  that  the  particles  are  approximately  spherical, 

m'"  = -coKC^/^lPeq  -  P\^  (15) 

where  a>=  .  This  result  can  be  inserted  directly 

into  the  cross-sectional  average  values  governing  Eqs.  (4)  and  (5), 
assuming  that  n'"  is  constant.  In  order  for  such  an  assumption  to 
hold,  the  following  conditions  must  be  satisfied: 

•  The  particles  are  evenly  distributed,  and  no  agglomeration  or 
break-up  into  smaller  particles  occur. 

•  Particles  do  not  accumulate  within  the  reactor,  nor  are  they 
completely  consumed  by  the  reaction. 

•  Particles  at  a  given  location  within  the  channel  are  subject  to 
the  same  local  pressure  and  temperature  conditions. 


The  energy  equation  for  the  slurry  side  and  load  side,  respectively, 
become: 

do 

=  -DaM"')!  -  a„)  +  (?)„(d„+,  -  d„)  +  (?)„_,  (d„-i  -  d„)  (17) 

do 

=  R„[(?)„(d„+1  -  dn)  +  (p„+,(e„_i  -  0„)]  (18) 

The  species  conservation  equation  becomes: 

^=-DaM';(l-a„)  (19) 

There  are  three  dimensionless  groups  that  arise:  a  Damkohler  num¬ 
ber  (Da),  (p,  and  heat  capacity  rate  ratio  {R).  The  Damkohler  number 
is  a  ratio  of  the  timescale  of  the  chemical  reaction  to  the  timescale 
for  mass  transfer  by  advection.  The  group  (p,  which  is  a  function  of 
axial  position,  is  a  dimensionless  heat  transfer  coefficient  that 
relates  the  magnitude  of  heat  transfer  by  wall  heat  flux  to  heat 
transfer  by  advection;  the  Reynolds  number  and  Prandtl  number 
are  evaluated  according  to  the  physical  properties  and  superficial 
velocity  of  the  slurry,  and  overall  thermal  resistance  is  used  in  place 
of  the  convection  coefficient  in  the  Nusselt  number. 

2.3.  Correlation  selection 

There  are  no  correlations  available  in  the  literature  developed 
specifically  for  a  reacting  multi-phase  flow  in  a  chevron  PHE.  The 
criteria  for  selecting  empirical  correlations  to  explore  the  thermo- 
hydraulic  performance  should  adhere  to  the  following  guidelines: 

1.  Conservative  estimates  are  provided  for  reactor  design. 

2.  The  correlations  are  applicable  to  a  wide  range  of  geometric 
parameters  for  PHEs. 


In  effect,  these  assumptions  yield  a  shrinking  core-type  model 


wherein  only  the  size  of  the  particles  change  as  the  reaction  pro¬ 
gresses,  with  the  reaction  rate  tending  toward  zero  as  the  surface 

area  available  for  reaction  also  tends  toward  zero. 

To  elucidate  the  important  parameters  influencing  the  solution 
of  the  system,  a  non-dimensional  analysis  is  performed.  Introduc¬ 
ing  the  following  dimensionless  variables  and  assuming  no  maldis- 

tribution  is  present: 

(16a) 

\Pls^p,ls  ) 

tln=^ 

^in 

(16b) 

II 

N 

(16c) 

Da  =  ^^|Pe,(V,„)-P„,J^ 

Cin  Uis 

(16d) 

Many  heat  exchanger  correlations  are  derived  from  heat  and 
momentum  transfer  data  (i.e.,/and  Nu)  averaged  over  the  entire 
channel  length.  However,  the  quantities  to  which  the  said  data  is 
correlated  will  vary  with  axial  position,  mostly  due  to  changes  in 
the  void  fraction  and  mass  quality  of  the  flow.  To  account  for  these 
variations,  it  is  assumed  that  these  correlations  may  be  applied 
at  the  prevailing  conditions  at  any  axial  location  z,  and  inte¬ 
grated  over  the  channel  length  -  in  essence,  a  quasi-equilibrium 
model. 

In  this  study,  a  semi-empirical  model  due  to  Martin  [19]  was  se¬ 
lected  for  the  single-phase  friction  factor,  since  it  correlates  data 
for  a  range  of  chevron  angles  fairly  well.  Martin’s  model  is  based 
on  observations  gathered  from  earlier  flow  visualization  studies 
of  the  flow  patterns  inside  PHEs  [20,21],  which  showed  that  two 
types  of  flows  exist  simultaneously:  crossing  substreams  that  fol¬ 
low  the  corrugations,  and  sinusoidal  wavy  flow  between  the  upper 
and  lower  plate  corrugations  as  the  fluid  travels  down  the  main 
flow  direction.  Martin’s  model  is  expressed  as  [19]: 

_ COS^ _  ^  1  -COSjg 

yj  ^0.18  tan  fi  +  0.36  sin  /?  +/o  cos  fi  i/3.8/] 


ff)  ^ 

De  RePr 

n  ^  (^P^p)ls 


M 


/// 

n 


Peq{Tn)-Pn  ] 

PeqiTld,in)-Patm\ 


2 


0 


Peq  ^  P 
Peq  ^  P 


(16e) 

(16f) 


where. 


/o 


64 

Re;, 


(1.81n(Reh)-1.5) 


-2 


Re^i  <  2000 
Re/i  >  2000 


/i 


i^  +  3.85 

39 

ReO.289 


ROh  <  2000 
Re^i  >  2000 


(21) 

(22) 


where  Re^  is  the  Reynolds  number  based  on  the  hydraulic  diameter, 
defined  as  =  4b/(^>AE,  with  being  the  area  enlargement  factor 
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given  by  the  ratio  of  the  actual  surface  area  to  the  projected  surface 
area  of  the  corrugated  chevron  plate. 

For  two-phase  flow,  no  reports  exist  on  models  or  correlations 
specific  to  the  reacting  multi-phase  flow  configurations  encoun¬ 
tered  in  this  study.  The  most  reasonable  example  in  the  literature 
akin  to  this  situation  is  that  of  air-water  vertical  upward  flow, 
which  a  limited  number  of  investigators  have  considered  [22- 
27].  In  the  reacting  slurry,  the  liquid  phase  and  the  gas  byproducts 
do  not  undergo  additional  phase  changes,  and  changes  in  the  flow 
quality  only  occur  through  chemical  reaction  rather  than  evapora¬ 
tion  or  condensation.  One  may  imagine  an  equivalent  air-water 
flow  in  which  gas  is  continuously  introduced  along  the  length  of 
the  channel;  the  implicit  assumption  in  this  analogy  is  that  the 
hydrodynamics  of  air-water  flow  is  similar  to  that  of  the  gas  evolv¬ 
ing  reacting  flow  considered  in  this  study. 

Tribbe  and  Miiller-Steinhagen  [23,24]  studied  air-water  two- 
phase  vertical  upflow  in  different  commercial  chevron  PHEs  having 
different  corrugation  depths  and  chevron  angles.  These  authors 
found  that  their  frictional  pressure  drop  data,  using  the  appropri¬ 
ate  single-phase  friction  factor,  could  be  well  correlated  by  the 
relation  [24]: 

a)TP  =  1.423  -  0.0724  lnX  +  1.031 /X  (23) 

where  X  is  the  Lockhart-Martinelli  parameter,  defined  as: 


X  = 


(24) 


Tribbe  and  Miiller-Steinhagen  [24]  claim  their  correlation  to  be 
valid  within  0.06  < X<  10.  From  the  experiments  conducted  in  part 
I  of  this  study,  the  frictional  pressure  drop  could  be  adequately 
described  by  this  correlation  when  a  correction  factor  of  3.0  was 
applied.  Since  the  reaction  rates  of  the  sodium  bicarbonate/acetic 
acid  system  are  much  greater  than  the  decomposition  of  AC  (since 
the  former  is  mixing-limited),  and  in  the  interest  of  obtaining  con¬ 
servative  estimates,  the  model  presented  here  does  not  apply  this 
correction  factor  to  Eq.  (23).  Likewise,  following  the  recommenda¬ 
tion  of  Kreissig  and  Miiller-Steinhagen  [22],  the  void  fraction  was 
estimated  using  Rouhani’s  correlation  [28]: 


a.  =  {[1  +  0.12(1  -  x„)](x„/p,,„  +  (1  -  Xn)/p,) 

Hg,n  ^ 

+1.18/«„,,„p°=)[g(7(p,  -  (25) 

where  the  total  mass  flux  and  quality,  respectively,  are  expressed 
as: 


Ktal  =  UisAPis  +  Cin,„(l  -  t]„)]  (26) 

=  (1  -  +  1  -  f7„)“’  (27) 

As  was  shown  in  part  I,  even  with  a  crude  estimate  of  the  reaction 
rate,  the  forgoing  equations  could  be  applied  using  local  values  for  r] 
and  pg  to  obtain  reasonable  (for  two-phase  flow,  at  least)  estimates 
of  the  total  pressure  drop. 

The  Nusselt  number  for  two-phase  flow,  lacking  any  correla¬ 
tions  specifically  for  heat  transfer  in  air-water  flow  in  chevron- 
PHEs,  requires  recourse  to  heat/momentum  transfer  analogies.  In 
particular,  Martin  [19],  Schliinder  [29],  and  Abu-Khader  [30]  have 
found  that  a  generalized  form  of  the  classical  Leveque  solution 
could  correlate  heat  transfer  data  obtained  from  chevron  PHEs  sur¬ 
prisingly  well.  Martin  modified  the  original  generalized  Leveque 
equation  to  fit  experimental  data  for  turbulent  flow  in  industrial 
PHEs  from  several  sources  to  obtain  the  following  expression: 

Nu,  =  0.U2PT^/\p/pJ^^\fRe^  sin(2/?)]“'"  (28) 


Martin  obtained  the  term  sin  {2p)  by  taking  D^/L  to  be  the  dis¬ 
tance  between  two  crossing  corrugations.  One  may,  plausibly,  ex¬ 
tend  Martin’s  adaptation  of  the  generalized  Leveque  solution  to 
the  case  of  two-phase  flow  by  recognizing  the  agitation  and  mixing 
caused  by  the  interaction  between  gas  bubbles  and  a  surrounding 
fluid  is  similar  to  that  of  turbulence  [31].  The  Leveque  solution  is 
an  asymptotic  solution  based  on  approximating  the  velocity  profile 
in  the  vicinity  of  a  solid  wall  as  linear,  a  more  involved  discussion 
of  which  may  be  found  in  Ref.  [32].  Since  the  only  change  here 
from  Leveque’s  original  analysis  is  the  inclusion  of  the  two-phase 
multiplier,  one  finds  that  NUh  =  NUhi^jpfRel).  Without  sufficient 
data  to  re-adjust  the  constants  in  Martin’s  adaptation  of  the  gener¬ 
alized  Leveque  equation,  we  opted  simply  to  insert  the  two-phase 
multiplier  into  the  bracketed  term  in  Eq.  (28): 

Nu,  =  0.122  [fisKRe^  sin(2/?)]  (29) 

The  exponent  0.374  and  leading  multiplier  0.122  are  expected  to  re¬ 
quire  adjustment  to  better  match  multi-phase  reacting  flow  in  a 
chevron  PHE.  As  long  as  the  thermal  boundary  layer  remains  thin, 
as  it  would  in  turbulent  flow  or  highly  agitated  two-phase  flow, 
one  can  expect  the  Leveque  correlation  to  be  an  appropriate  basis 
for  two-phase  flow  correlations  in  chevron  PHEs.  Eq.  (29)  implies 
that  the  added  presence  of  a  gas  phase,  which  will  yield  >  1, 
actually  enhances  heat  transfer.  Such  enhancement  was  observed 
by  Vlasogiannis  et  al.  [27]  in  their  study  of  heat  transfer  with  air- 
water  flow  in  a  vertical  chevron  PHE.  Based  on  their  experiments, 
Vlasogiannis  et  al.  reported  enhancements  in  the  values  of  heat 
transfer  coefficient  in  all  cases  where  air  was  injected  with  water 
at  the  inlet. 

2.4.  Numerical  procedure 

The  transport  equations  were  solved  using  the  finite-volume 
method,  with  each  channel  being  partitioned  into  equally  sized 
cells.  Eqs.  (7)  and  (17)-(19)  are  integrated  over  each  cell,  where 
upon  Gauss’  divergence  theorem  is  used  to  change  the  volume 
integrals  of  the  convective  and  diffusive  terms  into  fluxes  at  the 
cell  face.  1st  order  upwind  differencing  is  used  for  the  convective 
fluxes,  and  2nd  order  central  differencing  is  used  for  the  diffusive 
fluxes.  A  grid  resolution  of  100  cells  per  channel  was  found  to  be 
sufficient  to  establish  grid  independence. 

Addressing  the  nonlinear  source  terms  in  the  species  and 
energy  equations  require  an  iterative  approach.  The  energy  and 
species  conservation  equations  are  linearized  about  an  initial 
guess,  and  decoupled  to  form  two  separate  linear  systems,  which 
are  subsequently  solved  using  conventional  sparse  matrix  tech¬ 
niques.  The  pressure  field  is  updated  at  the  end  of  each  iteration. 
The  solver  code  was  written  in  Maple  13™,  and  the  matrix  inver¬ 
sion  was  performed  using  the  built  in  LinearSolve  function.  For 
the  simulations  presented  in  this  paper,  convergence  typically 
requires  50  iterations  and  3  min  of  wall-clock  time  for  execution. 


3.  Results  and  discussion 

The  numerical  model  described  above  was  used  to  design  and 
evaluate  the  performance  of  a  “prototype”  HEX  reactor.  The  reactor 
consists  of  7  channels,  3  of  which  carry  the  slurry,  and  4  of  which 
carry  the  thermal  load.  This  layout  is  shown  in  Fig.  3.  The  geomet¬ 
ric  parameters,  based  off  the  Alfa  Laval  M3  PHE  used  in  part  I  of  this 
study,  are  shown  in  Table  1.  The  intent  of  this  exercise  is  to  first 
determine  the  dominant  parameter  (or  parameters)  governing 
reactor  performance,  and  then  to  determine  if  the  operating  condi¬ 
tions  required  to  achieve  full  conversion  near  atmospheric  pres¬ 
sure  conditions  are  indeed  realistic. 
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Table  1 

Geometric  parameters  for  “prototype”  PHE. 

Parameter  Value 

L  0.48  m 

De  0.00468  m 

W  0.10  m 

^AE  1.1 

P  30° 


The  simulations  that  were  performed  in  pursuit  of  these  objec¬ 
tives  assumed  the  slurry  inlet  temperature  to  be  fixed  at  335  K, 
such  that  the  reaction  will  be  initiated  immediately.  The  mass 
loading  of  the  AC  was  held  at  8  kg/m^  while  systematically  varying 
the  slurry  flow  rate  and  inlet  pressure.  The  slurry-side  Reynolds 
number  ranged  from  5  to  16;  Abu  Kader  [30]  found  the  modified 
generalized  Leveque  equation  to  be  an  appropriate  Nusselt  number 
correlation  (in  single-phase  flow)  at  low  Re  for  the  chevron  angle  in 
Table  1,  but  no  data  exists  to  validate  Eq.  (29)  for  multi-phase  flow 
under  low  Re  conditions.  The  single-phase  friction  factor  correla¬ 
tions  in  Eqs.  (20)-(22)  are  also  applicable  to  laminar  flow,  and 
are  therefore  appropriate  for  the  flow  conditions  considered  in 
the  simulations.  All  thermophysical  properties  were  assumed  to 
remain  constant,  and  were  evaluated  at  the  slurry  inlet  tempera¬ 
ture.  This  particular  AC  mass  loading  was  chosen  based  on  a  target 
endothermic  heat  rejection  rate  of  2  kW,  assuming  near  total  con¬ 
version.  In  the  interest  of  obtaining  conservative  estimates  for  heat 
transfer,  the  viscosity  correction  term  in  Eq.  (29)  was  neglected 
(i.e.,  the  ratio  was  assumed  to  be  unity).  The  load  fluid  temperature 
is  held  constant  to  better  isolate  the  effects  of  slurry-side  convec¬ 
tion  and  reaction.  The  property  values  used  for  these  simulations 
are  shown  in  Table  2. 

The  volumetric  reaction  terms  in  Eqs.  (17)  and  (19)  require  an 
estimation  of  the  rate  constant,  k,  and  the  average  particle  size. 
As  an  approximate  first  estimate,  the  peak  initial  power  data 
reported  by  Schmidt  [7]  was  used,  along  with  the  assumption  that 
the  initial  value  of  the  average  particle  diameter  was  ~1  mm.  By 
performing  the  energy  balance  on  a  closed-system  -  representative 
of  the  reactor  used  by  Schmidt  -  and  using  the  reaction  rate  term 
from  Eq.  (15)  in  addition  to  assuming  the  slurry  temperature  is 
uniform,  the  instantaneous  power  rating  is  obtained  from 
w  =  AHAtotai^iP  -  Peq)^,  wheve  Afotai  is  the  total  particle  surface  area 
in  the  reaction.  Peg  is  calculated  at  the  fluid  temperature  using  Eq. 
(15),  and  P  is  assumed  to  be  the  set  point  pressure.  The  value  of  w 
was  estimated  from  the  peak  power  ratings  reported  by  Schmidt  at 
60  °C  with  a  set  point  of  400  torr  and  propylene  glycol  as  the  heat 
transfer  fluid.  The  reasoning  behind  this  assumption  was  that  the 
400  torr  setpoint  was  reached  quickly,  offering  a  more  accurate 
estimate  of  P;  and  the  maximum  power  occurs  early  in  the  exper¬ 
iment,  meaning  the  AC  is  closest  to  the  initial  reactor  temperature 
of  60  °C  and  thereby  allowing  more  accurate  estimation  of  Peg.  This 
calculation  yields  k:  =  2.4x  10“^°  m“^  kPa“^,  which  was  used  in 
the  simulations,  discussed  next. 


Table  2 

Thermophysical  properties  used  in  simulations. 


Property 

Value 

Units 

Source 

Opiate 

13.4 

W/mK 

[33] 

kis 

0.196 

W/mK 

[33] 

Cp,ls 

2700 

J/kgK 

[34] 

Ills 

8.4 

mPa  s 

[35] 

fig 

0.012 

mPa  s 

[33] 

Pis 

1036 

kgInP 

[35] 

Pp 

780 

kg/m" 

[6] 

Ah 

2000 

kj/kg 

[6] 

3.1.  Numerical  predictions 

One  of  the  major  concerns  for  open  flow  reactors  is  the  low 
residence  time  in  which  the  reaction  must  occur.  This  consideration 
places  a  limit  on  the  amount  of  AC  that  can  be  fed  into  the  reactor 
for  a  given  set  of  operating  conditions.  The  task,  then,  is  to  predict 
the  maximum  amount  (or  mass  concentration)  of  AC  that  can 
potentially  be  reacted.  There  are  three  major  competing  factors  that 
influence  the  tradeoff  between  conversion  efficiency  and  overall 
power  rating,  which  are:  residence  time/advection,  channel- 
to-channel  heat  transfer,  and  reaction  rate. 

Eig.  5  displays  conversion  as  a  function  of  load  side  temperature 
for  two  different  flow  rates  and  different  inlet  pressures.  Reducing 
the  inlet  pressure  from  110  kPa  to  102  kPa  did  not  have  a  signifi¬ 
cant  effect  on  the  load  temperature  at  which  the  AC  was  com¬ 
pletely  reacted,  whereas  increasing  the  slurry  flow  rate  by 
threefold  required  the  load  temperature  to  be  raised  by  ~10  K  to 
achieve  full  conversion.  The  conversion  efficiency  becomes  less 
sensitive  to  changes  in  the  inlet  pressure  at  higher  flow  rates. 
The  shape  of  the  conversion  curves  reveals  two  competing  effects 
on  the  reaction  rate;  surface  area,  and  equilibrium  drop.  The 
equilibrium  drop,  which  is  proportional  to  (Peg  -  Pf  increases  with 
axial  distance  due  to  decrease  in  local  pressure  from  hydrodynamic 
losses,  and  increase  in  equilibrium  pressure  as  the  slurry  tempera¬ 
ture  rises.  While  these  effects  tend  to  accelerate  the  reaction,  the 
decreasing  particle  surface  area  eventually  begins  to  dominate 
the  conversion  behavior,  causing  the  sensitivity  to  load  tempera¬ 
ture  to  severely  diminish  after  ~90%  conversion.  In  reality,  such 
asymptotic  behavior  is  unlikely  to  be  observed  since  particle  break 
up  and  changes  in  the  surface  morphology  would  have  an  increas¬ 
ingly  prominent  effect  on  the  total  solid  surface  area  per  unit 
volume  of  slurry. 

As  an  additional  figure  of  merit  for  assessing  reactor  perfor¬ 
mance,  it  is  insightful  to  consider  the  fraction  of  the  total  heat 
rejected  (Qioad)  in  the  form  of  latent  heat  (Qiatent)  as  a  supplement 
to  the  overall  power  rating.  Ideally,  the  proportion  of  heat  rejected 
to  the  thermo-chemical  heat  sink  should  be  as  high  as  possible.  In 
Eig.  6,  the  ratio  s  =  Qiatent/Oioad  is  plotted  with  respect  to  normalized 
load  temperature,  T*  =  (T/d  -  335  K)/(T|d  -  335  K),  where  is  the 
load  temperature  at  which  “full”  conversion  is  attained  (i.e., 
\  -rj  ^  0.95).  As  with  the  conversion  efficiency,  the  inlet  pressure 
has  little  effect  on  s  as  compared  to  the  residence  time.  Increasing 
the  flow  rate  tends  to  increase  the  total  power  rating,  but  at  the  ex¬ 
pense  of  8  since  higher  load  temperatures  are  required.  Thus,  the 


Fig.  5.  Comparison  of  conversion  efficiency  (h)  for  slurry  superficial  velocities  of 
0.01  m/s  (dashed  line)  and  0.03  m/s  (solid  line)  as  a  function  of  load  temperature  for 
different  inlet  pressures:  black  -  102  kPa,  green  -  105  kPa,  red  -  108  kPa,  and  blue  - 
110  kPa.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 
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(Tioad  "  Tsiurry)  !  (Tioad*  "^slurry) 


Fig.  6.  Fraction  of  total  heat  rejected  in  the  form  of  latent  heat  as  a  function  of 
normalized  load  temperature  for  different  inlet  pressures:  black  -  102  kPa,  green  - 
105  kPa,  red  -  108  kPa,  and  blue  -  110  kPa.  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 

increase  in  Nu  accompanying  lower  residence  times  is  outweighed 
by  the  higher  heat  capacity  rate. 

For  operating  conditions  close  to  atmospheric  pressure  condi¬ 
tions,  the  pressure  losses  are  insufficient  to  supply  the  equilibrium 
pressure  drop  needed  to  bridge  the  gulf  between  the  advection  and 
reaction  timescales,  as  evidenced  by  the  greater  load  temperatures 
required  at  higher  flow  rates.  It  may  be  inferred  from  this  fact  that 
the  global  reaction  rate  is  relatively  insensitive  to  the  axial  variation 
of  pressure  (at  least  at  the  flow  rates  considered  in  this  study),  and 
is  more  sensitive  to  the  axial  variation  in  the  slurry  temperature. 

From  the  standpoint  of  conversion  efficiency  and  heat  transfer,  the 
residence  time  is  the  dominant  parameter  for  operation  close  to 
atmospheric  pressure  conditions.  Increasing  the  residence  time  per¬ 
mits  a  greater  portion  of  the  total  heat  load  to  be  rejected  to  the  AC, 
but  the  overall  power  rating  is  diminished.  Additional  plates  may 
be  added  to  the  reactor  to  adjust  the  residence  time  to  the  appropriate 
level  if  one  desires  to  react  a  greater  quantity  of  AC  by  increasing  the 
slurry  flow  rate.  If  the  objective  is  to  maximize  the  proportion  of  the 
thermal  load  rejected  to  the  AC  for  a  fixed  PHE  design,  the  residence 
times  will  need  to  be  maximized  to  the  greatest  extent  possible. 

The  instances  where  e  >  1  correspond  to  conditions  when  the 
reaction  rate  outpaces  the  rate  of  channel-to-channel  heat  transfer. 
Hence,  the  reaction  consumes  sensible  heat  from  the  slurry  faster 
than  it  can  be  extracted  from  the  load.  Essentially,  the  AC  is  being 
wasted  in  this  condition  since  it  is  primarily  cooling  the  slurry 
rather  than  the  hot  stream  fluid,  as  is  evident  in  the  axial  temper¬ 
ature  profiles  in  Eigs.  7  and  8.  On  the  vertical  axes  of  Eigs.  7  and  8,  a 
new  dimensionless  temperature  is  defined  as  0  =  (T/s  -335  K)/ 

(Tid  -  335  K),  so  that  the  effect  of  the  endothermic  reaction  is 
readily  apparent.  The  solid  black  curve,  representing  the  case 
Tid  =  336  K,  can  be  seen  to  increase  before  gradually  decreasing 
when  the  inlet  pressure  is  110  kPa  (Pig.  7).  With  the  inlet  pressure 
at  102  kPa  (Pig.  8),  the  slurry  temperature  when  Tid  =  336  K  actu¬ 
ally  falls  below  its  inlet  value,  and  stabilizes  at  approximately 
334.8  K  before  slowly  increasing  again  due  to  the  increase  in  the 
heat  transfer  coefficient  as  the  reaction  progresses.  It  is  desirable 
to  have  s  as  close  to  unity  as  practically  possible,  since  heat  would 
be  extracted  from  the  load  fluid  at  the  same  rate  it  is  consumed  by 
the  endothermic  reaction. 

3.2.  HEX  reactor  sizing  calculations 

As  a  benchmark  for  HEX  reactor  performance,  the  data  obtained 
from  the  above  simulations  was  used  to  estimate  the  size  of  a 
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Fig.  7.  Dimensionless  axial  temperature  and  concentration  profiles  for  different 
load  temperatures  (black  -  336  K,  green  -  339  K,  red  -  342  K)  with  Uis  =  0.01  m/s, 
Qn  =  8  kg/m^,  and  =110  kPa.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  8.  Dimensionless  axial  temperature  and  concentration  profiles  for  different 
load  temperatures  (black  -  336  K,  green  -  339  K,  red  -  342  K)  with  Uis  =  0.01  m/s, 
Qn  =  8  kg/m^,  and  =102  kPa.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


reactor  and  the  required  operating  conditions  to  reject  2  kW  of 
heat.  Assuming  that  the  heat  rejection  per  channel  is  similar  to  that 
of  a  larger  reactor  (>20  plates).  Table  3  shows  the  minimum  reactor 
size  for  a  2  kW  heat  load  based  on  operating  conditions.  The  first 
four  entries  from  the  top  are  for  the  lowest  residence  time 
(Uis  =  0.01  m/s),  where  over  50%  of  the  total  heat  load  is  rejected 
to  the  AC;  with  30-35  plates,  the  Alfa  Laval  M3  would  occupy  a 
volume  of  ~30  L  with  a  footprint  of  865  cm^.  The  reactor  would 
be  roughly  the  same  size  as  a  typical  desktop  printer,  and  with 
an  AC  loading  of  8  kg/m^,  and  consume  about  0.7  g/s  of  solid  AC. 
With  the  slurry  temperature  at  335  K,  the  required  load  tempera¬ 
ture  is,  at  most,  342  K,  which  is  well  below  the  industry  standard 
safe  limit  for  electronics  cooling  (358  K)  [36].  The  plate  surface 
area  is  approximately  530  cm^  (the  actual  area  available  for  heat 
transfer  may  be  lower  since  the  corrugated  section  does  not  cover 
the  entire  plate),  from  which  the  heat  load  per  plate  q"  is  derived. 

Even  with  conservative  estimates  for  thermal  performance 
(neglecting  viscosity  correction  and  multiplier  correction  to  Eq. 
(23)),  the  size  and  operating  temperature  requirements  for  a 
small-scale  reactor  are  reasonable  for  operating  conditions  close 
to  atmospheric.  It  is  also  clear  that  the  tradeoff  between  residence 
time  and  total  power  rating  means  that  increasing  the  specific  heat 
duty  of  the  reactor  (heat  load  per  unit  system  volume)  causes 
depreciation  in  s.  In  order  to  strike  the  appropriate  balance  be¬ 
tween  s  and  the  overall  power  rating  will  depend  on  the  actual 
application  and  operating  constraints. 
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Table  3 

Reactor  size  (minimum  number  of  plates)  required  to  reject  2  kW  of  waste  heat;  AT  =  T^  -  335  K  is  the  load-slurry  temperature  difference  at  the  inlet. 


Slurry  flow  rate  (L/s) 

Qlatent  (W) 

Qtotal  (W) 

Qlatent  /Qtotal 

AT  (K) 

Pin  (kPa) 

Min.  plates 

q"  (W/cm^) 

0.078 

1114 

2017 

0.55 

7.4 

no 

30 

0.12 

0.083 

1185 

2076 

0.57 

6.9 

108 

32 

0.11 

0.083 

1185 

2029 

0.58 

6.9 

106 

32 

0.11 

0.09 

1293 

2040 

0.63 

5.8 

102 

15 

0.1 

0.049 

754 

2015 

0.37 

19 

110 

6 

0.58 

0.056 

862 

2275 

0.38 

18.7 

108 

7 

0.56 

0.056 

862 

2228 

0.39 

18.1 

105 

7 

0.55 

0.056 

862 

2168 

0.4 

17.5 

102 

7 

0.54 

4.  Conclusions 

In  order  to  provide  a  useful  engineering  analysis  tool  to  predict 
HEX  reactor  performance,  we  have  developed  a  reduced-order 
semi-empirical  model  since  the  analysis  of  multi-phase  flows  is 
still  a  heavily  empirical  practice  [11].  The  selection  of  empirical 
correlations  for  the  void  fraction  and  two-phase  multiplier  was 
guided  by  experimental  data  obtained  in  part  I.  As  an  initial  feasi¬ 
bility  study,  a  reduced-order  reactor  model  was  developed  to 
ascertain  the  critical  parameters  governing  the  performance.  The 
model  was  based  on  conventional  correlations  for  plate  heat 
exchangers.  The  numerical  results  show  that  the  conversion  effi¬ 
ciency  and  the  proportion  of  the  total  heat  transferred  consumed 
by  the  chemical  reaction  are  more  strongly  influenced  by  the  res¬ 
idence  time  than  the  operating  pressure.  Furthermore,  the  de¬ 
signed  reactor  meets  the  constraints  for  size  and  operating 
temperatures  while  achieving  the  required  thermal  load.  A  thermal 
load  of  2  kW  is  designed  to  be  handled  by  a  PHE  with  a  footprint  of 
only  30  L  in  volume  with  operating  conditions  close  to  atmospheric 
pressure. 
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